Strongly Correlated Superconductivity rising from a Pseudo-gap Metal 
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We solve by Dynamical Mean Field Theory a toy-model which has a phase diagram strikingly 
similar to that of high Tc superconductors: a bell-shaped superconducting region adjacent the Mott 
insulator and a normal phase that evolves from a conventional Fermi liquid to a pseudogapped 
semi-metal as the Mott transition is approached. Guided by the physics of the impurity model 
that is self-consistently solved within Dynamical Mean Field Theory, we introduce an analytical 
ansatz to model the dynamical behavior across the various phases which fits very accurately the 
numerical data. The ansatz is based on the assumption that the wave-function renormalization, 
that is very severe especially in the pseudogap phase close to the Mott transition, is perfectly 
canceled by the vertex corrections in the Cooper pairing channel. A remarkable outcome is that 
a superconducting state can develop even from a pseudogapped normal state, in which there are 
no low-energy quasiparticles. The overall physical scenario that emerges, although unraveled in a 
specific model and in an infinite-coordination Bethe lattice, can be interpreted in terms of so general 
arguments to suggest that it can be realized in other correlated systems. 

PACS numbers: 74.20.Mn, 71.27.-(-a, 71.30.-fh, 71.10.Hf 



I. INTRODUCTION 



How can high-temperature superconductivity emerge 
out of a pseudogap metal is one of the standing puz- 
zles posed by the cuprate superconductors. Indeed one 
would expect that a pseudogapped metal is not ideal for 
superconductivity mainly for two reasons: the pseudogap 
reduces the density of states at the Fermi level and, more 
worrisome, it is likely to cut off the BCS singularity in the 
Cooper channel. As a consequence, a very strong pair- 
ing should be required to turn such an unconventional 
metal into a high-Tc superconductor. Many alternative 
proposals have been put forward to reconcile the exis- 
tence of a pseudogap in the underdoped normal phase, 
which appears at a temperature T*, with the occurrence 
of superconductivity (SC) below a critical temperature Tc 
that may become significantly smaller than for deeply 
underdoped systems. 

The simplest explanation is to associate the opening 
of a pseudogap to the existence of preformed Cooper 
pairs well above Ti^-^^^i^. This is compatible with the 
low dimensionality and the high critical temperature of 
the cuprates, which cooperate to enhance phase fluctu- 
ations of the order parameter, leading to a wide region 
where pairs are already formed yet lacking phase coher- 
ence. This scenario is experimentally supported by the 
strong diamagnetic response observed in an extended re- 
gion above Tc^ . However, even though this point of view 
is certainly reasonable close to Tc, its application close 
to the pseudogap temperature scale and for small 
doping X is definitely more questionable. It is indeed 
well established that while increases monotonically as 
the doping a; — s- 0, both and the superfluid density 



vanish. Within the preformed-pair picture, this would 
correspond to very strong coupling leading to localized 
pairs in real space, which are hardly reconciled with the 
momentum-space nodal structure of the pseudogap ob- 
served by angle-resolved photoemissionil 

A different point of view interprets the pseudogap as 
due primarily to a competing ordered phas o^'^'^^d^d^d^ 
or arising from fluctuating competing orders among 
which c?-wave superconductivity prevails below Tc (see 
Ref. 14 and references therein). Although compatible 
with many experimental evidences, these proposals pose, 
in our opinion, several theoretical questions. For in- 
stance, even if we assume that the pseudogap phase is 
a fluctuating mixture of different orders, we are still left 
with the question about the nature of the underneath 
normal phase unstable to all the above competing or- 
ders. A common belief is that the antiferromagnetic 
ground state of the undoped parent compounds^ or an- 
other state very close in energj*i^ must be interpreted as 
the ancestor phases that naturally evolve upon doping 
into a novel state of matter— - a fluctuating mixture of 
pseudogap phases - rather than into a bona fide normal 
metal that, below Tc, turns superconducting. 

In this work we do not intend to enter these controver- 
sial issues in the context of cuprates. Rather, we want to 
unravel in all its facets a similar phenomenology - pseu- 
dogap normal phase which turns into a high-temperature 
superconductor - that we recently discovered by solving 
with Dynamical Mean Field Theory (DMFT)ii a two- 
orbital Hubbard model inspired by fullerene supercon- 
ductors.^* We think that providing an exhaustive anal- 
ysis of the pseudogap phenomenon in a strongly cor- 
related model that can be exactly solved, albeit in an 
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infinite-coordination lattice, may shed light on more re- 
alistic models for the cuprates, which are harder to deal 
with, both analytically and numerically. We are going 
to show that in this two-orbital model, akin to models 
for cuprates, superconductivity is the low-tcmperature 
winner among competing phases. In our case, two of 
these competing phases are homogeneous and symmetry- 
invariant: a conventional Fermi-liquid metal and an in- 
trinsic pseudogap phase, i.e., a single quantum phase 
with zero entropy at T = 0. The connection between 
the superconducting and the intrinsic pseudogap phases 
is intriguing, as the latter proves to be a fertile ground 
for superconductivity. The analysis of the self-energy in 
the superconducting phase reveals that the onset of su- 
perconductivity gives rise to much more coherent excita- 
tions than in the normal pseudogap state, where there are 
no quasiparticles and the self-energy diverges at low fre- 
quency. The way in which this singularity is regularized 
in the superconducting phase bears strong similarities 
with the effect of non magnetic impurities in an s-wave 
superconductor, where the processes leading to anomalies 
in the normal phase do not cut-off the superconducting 
instability, or, in other words, they are non-pairbreaking. 
Despite the regularization of the low-frequency anoma- 
lies, where the superconducting gap develops, our corre- 
lated superconductor still presents the pseudogap energy 
scale at high-frequency, and these two scales have oppo- 
site behavior as the doping goes to zero. 

Although we are going to review quite in detail the 
model and its properties in the following sections, we 
think is worth anticipating here some of the main results. 
In Fig. [1] we sketch the phase diagram of the model. 
On the X-axis we plot a parameter that controls the dis- 
tance from a Mott insulating phase like doping or pres- 
sure. The acronyms FL, NFL, SC and PG stand for nor- 
mal Fermi-liquid, non-Fermi-liquid, superconducting and 
normal pseudogap phases, respectively. The temperature 
scale r_ (two branches) identifies the crossovers from the 
FL and PG phases into the almost critical NFL region. 
The latter is not pseudogapped and it is characterized 
by incoherent single-particle excitations with dispersion 
and inverse-lifetime controlled by a single scale T^. The 
Mott transition occurs when T-|_ = r_ on the pseudogap 
side of the superconducting region. 

The phase diagram Fig. [1] bears striking similarities to 
that of cuprates, identifying the T- branch on the "un- 
derdoped" side with . From a purely theoretical point 
of view, the diagram closely resembles the quantum- 
critical-point (QCP) scenario proposed in many con- 
texts, including cuprates^i^iii^ and heavy-fermion com- 
poundsiiii^ii^, as well as the phase diagram of the t — J 
model for cuprate a^"^i^^'^° within the Resonating Valence 
Bond (RVB) frameworkfi^ We postpone a critical com- 
parison with those to the ending section of the paper. 

The paper is organized as follows. In Section II we in- 
troduce our model and approach. Sec. Ill is devoted to 
the DMFT phase diagram, while Sec. IV summarizes the 
information that we can obtain from the impurity model 
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FIG. 1: (Color online) Sketch of the phase diagram of the 
model Eq. ([T}. a; is a parameter that measures the deviation 
from the Mott insulator. FL, NFL, SC and PG stand for 
normal Fermi-liquid, non-Fermi-liquid, superconducting and 
normal pseudogap phases, respectively. T_ and T+ are cross- 
over lines, see the text. 

which corresponds to our lattice model through DMFT. 
Sees. V and VI introduce modelizations of the self-energy 
in the normal and superconducting phases, respectively. 
Finally, Sec. VII is dedicated to conclusions with par- 
ticular attention to the relation between our results and 
theoretical proposals for the cuprates. 



II. THE MODEL 

The model that we study is a two-orbital Hubbard 
model with the Hamiltonian 

2 

a— 1 G 

i 

-2JY,{tI^+TI^, (1) 

i 

where c| and Cj creates and annihilates, respectively, 
one electron at site i in orbital a = 1, 2 with spin cr, 
ni ~ ni^\+ riiji is the on-site occupation number, being 

n^^a = Y.aA,aa%aa^ ^ud finally, 

^i," — 2 E/ $Z '^itO." ^ab bo- (2) 
a,&=l (J 

are orbital pseudo-spin- 1/2 operators, r" (a = x,y,z) 
being the Pauli matrices. We assume hereafter that the 
exchange constant J is positive, hence favors low-spin 
atomic configurations. This model was introduced in 
Ref. 18 to mimic an. eiS)E Jahn- Teller coupling to a local 
doubly-degenerate phonon mode that prevails over the 
conventional Coulomb exchange. In this case, the Jahn- 
Teller coupling leads effectively to inverted Hund's rules. 
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on provision that the phonon frequency is high enough 
to neglect retardation effects. The original purpose was 
to study a simplified model that shared the same physics 
of alkali-doped fuUerene superconductors, where pairing 
is mediated by eight fivefold-degenerate local vibrational 
modes, t(^H Jahn- Teller coupled to the threefold degen- 
erate LUMO of Cgo'^'^'^ Apart from a constant term, 
the Hamiltonian ([T]) can be alternatively written as 

2 

ij a—1 a 
i a—1 

+ E "^J- ■ + Vni4 ni,2, (3) 

i 

with = 4J and F = C/ + J, that also describes two 
Hubbard models labeled by the orbital index a = 1,2, 
coupled by an antiferromagnetic exchange Jj_ and by a 
strong repulsion V . 

The interaction term proportional to J in ([1]) is easily 
shown to generate an attraction that leads to an order 
parameter of s-wave symmetry associated to the operator 

The bare scattering amplitude in this channel \s A — 
— 2J, and it would induce a superconducting instability 
in the absence of Coulomb repulsion. The introduction 
of U weakens the attraction, leading to A = — 2J -I- U . 
Therefore, at least in weak coupling, l],J (W^ being 

the non- interacting bandwidth), one expects the model 
to describe a BCS superconductor for U < 2 J, and a nor- 
mal metal for larger repulsion. In principle, if the non- 
interacting Fermi surface accidentally has nesting, other 
weak-coupling instabilities like, e.g., magnetism can oc- 
cur. In our calculation we did not consider such commen- 
surate phases, because they are related to specific aspects 
of the lattice, and we want to focus on more basic and 
general properties. We note that the competition be- 
tween attraction and repulsion, namely between J± and 
1^ in ([3]) , offers the opportunity to investigate how s-wave 
superconductivity can emerge at all in spite of strong 
correlations. This issue is not commonly touched in the 
context of cuprates, where the emphasis is mostly put on 
the appearance of d-wave superconductivity (whose order 
parameter lives on the bonds) in the presence of purely 
on-site repulsion, while longer-range contributions to the 
interaction are often ignored. 

Let us move now to the opposite limit of strong re- 
pulsion, U ^ W, and consider the half-filled case. In 
this limit the model is a Mott insulator, each site being 
occupied by two electrons that can not move coherently. 
When J — the model maps onto an SU(4) Heisenberg 
model. If the hopping is restricted to nearest neighbor, 
the ground state of this model is dimerized in one di- 
mension^^, while in a two dimensional square lattice is 



still unclear whether is spin-liquid or Neel orderedi ^^i^ 
When the attraction is switched on and J ^ W'^/U > 
the attractive interaction prevails and two electrons on 
each site i lock into the singlet state 

^ (4iTcl2i+cl2TcIu) |0>, (5) 

stable to the weak inter-site superexchange. In this case, 
the Mott insulator is non-magnetic and translationally 
invariant - a local version of a valence-bond crystal - 
regardless the structure of the hopping matrix element 
and the dimensionality and topology of the lattice. A 
sufficient degree of frustration, either due to the lat- 
tice topology or to the hopping amplitudes t^'s, dis- 
favors a magnetic state and makes it possible for this 
non-magnetic phase to survive decreasing U/W down to 
the Mott transition, which is in turn pushed to a finite 
Uc when frustration eliminates nesting. When this hap- 
pens, one should naively expect, increasing U/W at fixed 
J/W <C 1, first a very narrow BCS superconducting re- 
gion for < U < 2 J, followed by a normal metal that 
eventually gives way to a non-magnetic Mott insulator 
when U > Uc ^ W. Seemingly, doping the non-magnetic 
Mott insulator at U > Uc ^ J should bring to a nor- 
mal metal that, as doping increases, gets less and less 
correlated. 

This nai've expectation, based essentially on the value 
of the bare scattering amplitude A — —2J -I- J7 in the 
singlet Cooper-channel Eq. ([4]), turns out to be wrong at 
least in the two cases that have been so far considered: (i) 
an infinite-coordination Bethe lattice^*, which is exactly 
solved by DMFT; (ii) a one dimensional chain3^. We will 
briefly mention the latter in the last Section, while in 
what follows we concentrate mostly on the DMFT results 
for the Bethe lattice. In the following we first recall the 
basic ideas behind DMFT^ - and our implementation. 

DMFT extends the idea of classical mean-field to the 
quantum domain: a lattice model is approximately solved 
by solving the quantum problem of a single site subject to 
a "dynamical Weiss field" that describes the action of the 
rest of the lattice sites on the given site (assumed to be 
equivalent to any other). As in classical mean-field, the 
mapping is exact only in infinite coordination-lattices. 
The effective action of the local degrees of freedom reads 

Seff = f drdr' 4_(r) G,' (r - r')f , cop.'{r') + 
Jo 

+ 5o[coa<T,cJ„„], (6) 

where Qq^ is the Weiss field, a and P are orbital indices 
and a and cr' are spin indices. Sq is the local part of 
the action and contains all local interaction terms of the 
lattice Hamiltonian. The mean-field scheme is closed by 
imposing between the Weiss field and the local Green's 
function G (r — r'), computed by the action a self- 
consistency relation that contains the information about 
the original lattice model through the non-interacting 
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density of states (DOS), For a Bethe lattice with nearest 
neighbor hopping and bandwidth W, which we consider 
hereafter, the self-consistency reads 



Gq ^[i^nTa^, = 



E 



,af3 



W n 



(7) 



where E is the matrix of the single-particle terms 
of the on-site Haniiltonian (chemical potential, mag- 
netic field, hybridization) and G{iLUn) is the Fourier 
transform of the local Green's function G'(t)^^, — 
- {TTCoaa{T)clf^„, {Q)) s^f f ■ Eqs. © and can be viewed 
as a set of two coupled equations for G and Gq^ ■ In prac- 
tice, one needs to solve ^ for a given choice of Gq^ and 
obtaining G. Using ^ one finds a new value of Gq^- The 
procedure is iterated until convergence is achieved. It is 
evident that the computation of G is the hard part of the 
calculation. An important observation is that the effec- 
tive local theory can be represented as an impurity model 
whose hybridization function coincides with the dynami- 
cal Weiss field. In practice, a DMFT calculation amounts 
to solve an Anderson impurity model iteratively in order 
to self-consistently determine its hybridization function. 

The solution of the impurity model requires either an 
approximate numerical method, or an "exact" numerical 
approach. In this work we use exact diagonalizaton at 
T — 02^, which amounts to approximate the continuous 
bath of the impurity model by a discrete set of energy lev- 
els hybridized with the impurity. The method has been 
shown to converge exponentially as function of the num- 
ber Nb of bath levels. For example, Ni, — 5 already gives 
quite reliable results for the phase diagram and thermo- 
dynamic observables of a single band Hubbard model. In 
this work we will use 6 bath levels for each orbital, which 
gives Nb = 12 in total. An important aspect of the exact 
diagonalization method is the way in which the continu- 
ous bath is approximated. This is implemented by mini- 
mizing a suitably chosen distance between two functions. 
The distance is typically computed on the imaginary axis 
on a Matsubara grid corresponding to an effective tem- 
perature P, which is not to be confused with the physical 
temperature, always set to T = 0. In this work we typi- 
cally use (3 — AOO/W, and define a distance that weights 
more low frequencies'^'^. Specifically, we minimize 



E 



discrete \ 



UJ„ 



(8) 



{Go^)discrete bciug the discrctc version of the Weiss field, 
with the sum extending up to a maximum frequency of 
order AU. 

Finally, in order to describe the superconducting 
phase, the bath includes superconducting terms lead- 
ing to a Weiss field with anomalous components. Con- 
sequently the impurity Green's function has also an 
anomalous term Fa/sir) = —{TrCa'\{T)cf3i{0)), which is 
used together with G to build the matrix G(itUn) in 
the Nambu-Gor'kov formalism. Seemingly one define 
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FIG. 2: (Color online) DMFT phase diagram of the model 
Eq. ID) as function of U/W and doping 5. SCS stands for 
strong-correlated superconductivity. 



a matrix Weiss field with diagonal, tj"(ici;„)^^, and off- 
diagonal, T^{iujn)~^, components. The self-consistency 
condition ([7]) can be rewritten in the same formalism 
straightforwardly. 



III. THE PHASE DIAGRAM IN A BETHE 
LATTICE 

In Fig. [2] we show the DMFT phase diagram of model 
^ at fixed J = 0.05 VF in a Bethe lattice as func- 
tion of doping and U/W around the Mott transition, 
Uc — 0.87 WJ^. The first remarkable thing to note is 
the appearance of superconductivity, denoted as SCS (for 
Strongly Correlated Superconductivity, see below) in the 
figure, just around the Mott transition. The symmetry 
of the order parameter is that of Eq. ([4]). Superconduc- 
tivity even extends for U > Uc within a finite doping 
interval and is preceded by a pseudogap metal and fol- 
lowed, at larger doping, by a more conventional normal 
metal, see Fig. [2] We emphasize that the discreteness 
of the spectra obtained in ED does not allow to unam- 
biguously identify the pseudogap region, yet the results 
clearly show the evolution from one kind of normal state 
into the other, as already shown in Ref. [Tl . 

The strength of pairing is also surprising. In Fig. [3] we 
draw the values of the superconducting gap A at half- 
fiUing as function of U /W, panel (a), and at U = 0.92 W 
as function of doping, panel (b). In panel (a) we also 
show on a smaller scale the same quantity in the BCS- 
like region for small U < 2J. The latter is expo- 
nentially small, in agreement with the BCS estimate 
A ~ 0.5W exp {-ttW/SJ) ~ 2 x lO'^W, almost two 
orders of magnitude smaller that the values attained in 
the superconducting region around the Mott insulator. 

This superconducting phase that re-emerges with 
strengthened pairing just before the Mott transition was 
named in Ref. l23| Strongly Correlated Superconductiv- 
ity (SCS) to emphasize its peculiar properties with re- 
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FIG. 3: (Color online) Superconducting gap at half-filling as FIG. 4: (Color online) (a) Drude weight at half-filling as func- 

function of U/W, panel (a), and for U — 0.92 W as function tion of U/W. In the SCS region is also shown the Drude 

of doping away from half-filling, panel(b). The inset in panel weight of the metastable normal solution, zoomed in the in- 

(a) shows the gap at weak U on an enlarged scale. set; (b) Drude weight at [/ = 0.92 W as function of doping. 



spect to conventional BCS superconductors. Indeed, be- 
sides the large value of the gap in spite of the tiny at- 
traction J 0.05 W, other features characterize SCS. 
In Fig. IHa) we compare the values at half-fiUing of the 
Drude weight (zero-frequency contribution to the optical 
conductivity o'iuj)) in the SCS phase and in what could be 
regarded as the normal phase, namely the metastable so- 
lution obtained within DMFT by preventing gauge sym- 
metry breaking. We note that the onset of superconduc- 
tivity is accompanied by an increase of Drude weight, 
unlike what happens in a BCS superconductor. Remark- 
ably, while the SCS Drude weight vanishes only at the 
Mott transition, the weight of the normal solution van- 
ishes for a smaller U when the zero-frequency spectral 
weight goes to zero, suggesting the opening of the pseu- 
dogap. This semiconducting normal-phase is metastable 
at half-filling, where the stable zero-temperature phase 
is superconducting, but it is stabilized away from half- 
filling as shown in the phase diagram of Fig. [2l Here, 
doping seems to have the same effect as in semiconduc- 
tors, leading to a Drude weight linearly increasing with 
the number of holes, reported in panel (b) of Fig. 01 The 
intrusion of SCS leads to a more pronounced increase 
of Drude weight that smoothly connects to the normal 
metal appearing for larger doping. 



perturbation theory, the effect of interaction on the low- 
energy single-particle properties can be absorbed into the 
so-called wave-function renormalization or quasiparticle 
residue Z defined by 



Z \ diujn 



(9) 



where S(ia;„, kp)\s the single-particle self-energy in Mat- 
subara frequencies at the Fermi momentum, and by the 
reduction of the quasi-particle bandwidth W Wi, < 
W, 



W 



Z 



4 dk 



k—kp 



(10) 



where Vp is the bare Fermi velocity. Z can be regarded as 
the component of the total single-particle spectral weight 
that is carried by coherent quasiparticle excitations. In 
general Z < W^jW , the equality holding only in infinite- 
coordination lattices where the self-energy is momentum 
independent, E(iu;„,/c) = S](ici;„).'^l 

Within Landau theory, considering a generic scattering 
channel with bare amplitude A, the renormalized value 
can be written as 



A. Fermi liquid description 

Since the stable metallic phases at small U jW at half- 
filling, or at large doping away from half-filling, do not 
show any particular anomaly, a Landau Fermi-liquid sce- 
nario is expected to be applicable to understand how su- 
perconductivity can emerge in spite of the fact that the 
bare scattering amplitude in the Cooper channel is re- 
pulsive. The Fermi-liquid behavior is indeed confirmed 
by the regular self-energies calculated by DMFT. Within 



(11) 



where Fa includes the so-called vertex corrections. Ap- 
proaching an interaction-driven Mott transition, [/ — > 
t/c, the quasiparticle residue vanishes Z oc J7c — C/ ^ 0, 
but the behavior of A* in different channels can be to- 
tally different due to a different relevance of vertex cor- 
rections. Physical intuition suggests that the proximity 
to a Mott transition affects primarily charge fluctuations, 
that are severely suppressed, but it does not equally influ- 
ence the way in which the charge is distributed between 
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different spin and orbital states. Indeed, the localization 
of the charge leads to the formation of local moments 
and reflects an enhancement of the spin and orbital re- 
sponses. This suggests that, while vertex corrections in 
the spin- and orbital-density channels can compensate 
the vanishing Z, the same cancellation does not take 
place in the charge-density channel. Following this ob- 
servation, we argued in Ref. ,23 that the on-site repulsion 
U and exchange J undergo different renormalization as 
the Mott insulator is approached. In particular, since 
the exchange term J only controls the multiplet split- 
ting at fixed charge, it is not affected by the approach 
to the Mott transition, hence J ^ J* ~ J. On the con- 
trary, the residual quasi-particle repulsion is substantially 
weakened near a Mott transition, since most of correla- 
tion effects are already built into the small Z. On the 
basis of the DMFT behavior of the double occupancy in 
the single band Hubbard model, it was speculated in 
Ref. [23 that U —^U^~ZU. This assumption would 
result in the ansatz for the quasiparticle scattering am- 
plitude in the Cooper channel (jH) 



0.15 



A = - 2 J ^ A* = - 2 J* ~ Z [/ - 2 J. 



(12) 



Should Eq. ^2|) be correct, it would imply that A^,, which 
for small U > 2 J is repulsive in agreement with perturba- 
tion theory, must necessarily turn attractive sufSciently 
close to the Mott transition, where Z goes to zero. In 
Fig. [5] we show given by Eq. (|12p with Z extracted 
according to Eq. ^ from the normal-state DMFT so- 
lution. We note that this estimate for A, changes from 
repulsive to attractive very close to the point where a 
stable superconducting solution is found, supporting the 
validity of Moreover, the expression (fT^ provides 

an explanation for the large strength of pairing in the 
SCS phase in comparison with the BCS state. In fact, as 
Z ^ 0, a, regime in which A^ ~ — = —ZW is even- 
tually reached. Here the quasiparticles feel an attraction 
of the same order of their effective bandwidth. This in- 
termediate regime, bridging between the weak-coupling 
BCS limit and the strong-coupling Bose regime, has been 
shown to be the optimal situation for superconductivity 
in purely attractive models'^^'^'^. 

According to the above Fermi- liquid arguments, su- 
perconductivity re-emerges before the Mott transition 
because our Hamiltonian contains a pairing mechanism 
whose only role is to favor the singlet configuration ([5]) 
whenever two electrons sit on the same site. For this 
reason pairing is strengthened rather than weakened as 
Mott localization is approached. A similar phenomenon 
is not expected to occur if pairing is mediated by a 
charge-charge attraction, that obviously conflicts with 
the Coulomb repulsion. Notice that in the present con- 
text the origin of J (purely electronic generated, e.g., by 
superexchange, or driven by coupling with Jahn- Teller 
phonons as in fuUerenes-l) is not important; the only 
relevant point being that it does not compete with U. 
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FIG. 5: (Color online) Ansatz for the scattering amplitude 
A, = Z U — 2J versus U/W using Z extracted from a nor- 
mal DMFT solution. Also shown is the superconducting gap 
multiplied by a factor 20 for graphic convenience. 



IV. INSIGHTS FROM THE IMPURITY MODEL 

Although Fermi-liquid theory can be rather safely in- 
voked to explain the emergence of superconductivity from 
the correlated metal, it has to be eventually abandoned 
in the region of coupling which precedes the Mott transi- 
tion, where, as we discussed, superconductivity emerges 
out of a pseudogap normal phase. This demands for 
an alternative description able to account for both the 
Fermi- liquid and non- Fermi-liquid regimes. 

As we introduced in Sec. HH DMFT establishes a cor- 
respondence between a lattice model and an impurity 
model. The equivalence between the two models is en- 
forced by a self-consistency condition which contains the 
information about the original lattice. The role of the 
DMFT self-consistency is absolutely non trivial: for in- 
stance it makes the impurity Kondo temperature Tk van- 
ish at a flnite value of the interaction which signals the 
Mott transition in the lattice model, since Tk coincides 
with the renormalized bandwidth = Z W. How- 
ever, important insights can be gained by the analysis 
of the impurity model alone, without imposing the self- 
consistency. 

The Hamiltonian of the impurity model corresponding 
to H]) is 



k a— 1 (7 



u 
'2" 



(13) 



where cj^. and c^. are auxiliary fermionic degrees of 
freedom introduced to describe the Weiss fleld, while d\^, 
rfocri ^^"i T^(y)(z) ^re the same operators defined above, 
but they are specialized to the impurity site. The physics 
of the model is controlled by three energy scales: U, J 
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and the so-called hybridization width T defined through 
r^nJ2\Vw\' Sie^). (14) 

k 

The Hamiltonian (|13p has been recently studied in 
Refs. I33I35II36I by means of Numerical Renormalization 
Group (NRG) for a constant hybridization functioni^ii^ 
Analyzing the evolution of the physics as a function of 
U at fixed T and J <C F, a quantum phase transition 
takes place at U = C/,, an interaction value that corre- 
sponds to a Kondo temperature ~ J. For U < U^., 
i.e., Tk > T^, perfect Kondo screening takes place. On 
the contrary, when U > U^, Tk < Tj^, the impurity locks 
into the singlet configuration 

(4^4^ + 4^dl^) io), (15) 

which does not require any Kondo screening. The crit- 
ical point that separates the two phases is similar to 
that found in the two impurity Kondo model^i^. More- 
over the critical point exists also away from particle-hole 
symmetry, and, remarkably, it can be accessed even for 
U > Uf: changing the density away from half-filling.'^^ 
The single-particle spectral function displays an inter- 
esting evolution across the transitioni^^ In the Kondo 
screened phase, the low-frequency DOS is characterized 
by a broad resonance around the Fermi level on top of 
which a narrow Kondo peak develops. The latter shrinks 
continuously as the critical point is approached, while 
the broad resonance stays unaffected as well as the high- 
energy Hubbard bands surrounding the low-energy fea- 
tures. At the critical point, the Kondo peak disappears 
leaving behind only the broad resonance. In the un- 
screened phase, a narrow pseudogap appears inside the 
broad resonance, whose width grows moving away from 
the critical point. This pseudogap is gradually filled away 
from particle-hole symmetry, although, as mentioned, 
the critical point still exists. This behavior has been 
parametrized^ by the following ansatz for the low-energy 
DOS at particle-hole symmetry: 

pW = ^(^±^)^ (16) 

where the plus and minus signs refer to the screened 
and unscreened phases, respectively. The energy scale 
T_|_ ^ max{Tx^ J) measures the width of the broad reso- 
nance, while T- cx |?7 — J7*p controls the deviation from 
the fixed point and sets the width of the Kondo peak 
in the screened phase and of the pseudogap in the un- 
screened one. This model DOS defines a self-energy that 
was shown to fit perfectly well the numerical data for the 
impurity model.- » In particular, the self-energy at small 
Matsubara frequencies is Fermi-liquid-like in the Kondo 
screened phase, T,(iujn) oc —iujn, is finite and imaginary 
at the critical point, T,{iuJn) = — iF, and diverges in the 
unscreened phase, Ti{iuj„) cx 1/iuJn- 



The critical point is unstable in several symmetry 
breaking channelsS^^.^ including the magnetic channel 

S1-S2, (17) 

the Cooper channel Eq. ^ , and the hybridization chan- 
nels 

5] dld,^, J2 did,,, (18) 

that break the 0{2) orbital symmetry. All these chan- 
nels are degenerate at the critical point, where the model 
has an enlarged SO {7) symmetry.— When the impu- 
rity model contains explicitly a symmetry-breaking term 
that couples to one of the unstable channels, the crit- 
ical point is washed out and the phase transition be- 
tween the Kondo and the local singlet phase turns into 
a sharp cross-over , '^^1'^^ In the meantime, the relevant 
perturbation cuts off the non-Fermi-liquid singularities 
of the self-energy both at the critical point and in the 
unscreened phase, so that, at sufficiently low frequency, 
E(ia;,i) (X —iujn is always recovered;^ 

Since the Kondo temperature Tk has to vanish as 
the Mott transition is approached, the effective impurity 
model must encounter the critical point, Tk — J, before 
the Mott point is reached. In Ref. [s^ it has been spec- 
ulated that, once a full DMFT calculation is carried on, 
the instabilities associated to the impurity critical point 
lead because of the self-consistency condition to a spon- 
taneous symmetry breaking in a whole region around the 
critical point along one of the instability channels. These 
are the particle- hole channels ^T7\ and (fT5|) , which corre- 
spond to magnetic and orbital ordering, respectively, and 
the particle-particle channel Q that implies supercon- 
ductivity. This prediction is perfectly compatible with 
the DMFT phase diagram. Fig. [21 We note that particle- 
hole instabilities usually require nesting or other band- 
structure singularities that only accidentally occur, while 
the Gooper singularity is more ubiquitous. For this rea- 
son, in Ref. ilSj we only searched for a superconducting 
instability, even though the Bethe lattice with nearest 
neighbor hopping has nesting at half-filling. 



V. MODELING THE DYNAMICS IN THE 
NORMAL PHASES 

The analysis of the previous section shows that the 
critical point of the impurity model is expected to be al- 
ways pre-empted by broken symmetry phases in the lat- 
tice model treated within DMFT. Nonetheless, it is plau- 
sible that the critical point correspond to a metastable 
phase, just like a normal metal is a metastable phase in 
the presence of attraction. Following this idea, we have 
solved the model ([T]) forcing the DMFT self-consistency 
not to break any symmetry. The behavior of the self- 
energy that we obtain, see Fig. [HI closely reminds that 
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FIG. 6: (Color online) iti;„ — E(jaj„) versus Wn > for different 
values of U/W, as obtained by DMFT preventing supercon- 
ductivity. 



of the impurity model with constant bath (i.e. with- 
out self-consistency) that we just discussed. This con- 
firms that the impurity critical point corresponds to a 
metastable phase in the lattice model and suggests that 
the same parametrization Eq. (fT6|) may work in the lat- 
tice as well. Therefore we have assumed for the low fre- 
quency local Green's function the following expression, 
valid in a Bethe lattice (all energies will be expressed in 
units of the bandwidth W — 1), 



Go 



t1 



±50 



lUJri 



(19) 



where the ± sign, r+ and r_ have the same meaning as 
in Eq. (HU, while 



LOn - sign(w„) A / ^2 -I- 



(20) 



is the non-interacting local Green's function. The ansatz 
()19p for the Green's function corresponds to an ansatz 
for the self-energy which, through the DMFT self- 
consistency equation, becomes 

E(iLJ„) = iuJn - -^GiiLUn) - ^/(zcJn)"^ (21) 

16 

In the Fermi liquid region, where the plus sign has to be 
used in p^ and r_ ^ 0, the low-frequency self-energy is 



S(icj„) 



1 



2 \T+ 

,2/1 I 



1 

sign(c^„), (22) 



corresponding to a regular Fermi liquid behavior with 
quasi-particle residue 



t~Tt1 



(23) 



At the critical point (T_ = 0), or for frequencies T_ ^ 
<C T+, the self-energy becomes 



. 3 . . 5w„ 



sign{ujn). 



(24) 



E(iw„) has a finite and sizeable imaginary part, implying 
a non- Fermi liquid behavior. The deviations from con- 
ventional Fermi liquid behavior become even more pro- 
nounced in the unscreened phase, where T_ ^ and the 
minus sign has to be used in (fT9|) . Here we obtain 



4w„ T+ - T 
i T+ + T. 



4 T4. 



sign(w„) - i 



(25) 



The divergence of E (jw„) for small a;„ leads to a pseudo- 
gap in the DOS, whose behavior at small energy is 



1 



1 



(26) 



We can check the validity of our ansatz by simply com- 
paring it with the actual DMFT results obtained using 
exact diagonalization. Fitting these data with p9)l we 
find a very good agreement, and we can extract the pa- 
rameters Tj^ and r_ . The behavior of the best-fit values 
of T+ and is drawn in Fig. [7|and supports our ansatz, 
and consequently the relevance of the impurity critical 
point for the lattice model. In the same figure we report 
the superconducting gap obtained in DMFT allowing for 
gauge-symmetry breaking. 
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FIG. 7: (Color online) Fitting parameters r+ and T_ of the 
DMFT normal self-energy with the expression Eq. (ansatz- 
DMFT-Sigma) using for the impurity Green's function 
Eq. p9p . We also draw the DMFT superconducting gap when 
gauge symmetry breaking is allowed. 

Two things are worth noting in this figure. First, the 
vanishing of T-^ at C/ ~ 0.82 W , namely the location of 
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the metastable critical point, corresponds to the maxi- 
mum of the superconducting gap once gauge symmetry 
breaking is allowed, in accordance with the prediction 
based on the analysis of the impurity modeL-^ At the 
metastable critical point, the only remaining energy scale 
is T-f that therefore determines the maximum value of 
the superconducting gap A. Secondly, the Mott transi- 
tion occurs within the pseudogap region when T_ = T-|_. 
In other words, the pseudogap gradually widens as U in- 
creases from the T_ = point and, at the same time, 
the total spectral weight of the low-energy part, given 
by (r_|. — TL) according to Eq. (|19p . diminishes until it 
disappears, when r_ = r+. 

Although we solved DMFT at zero temperature, we 
believe that the T — energy scales will at least ap- 
proximately determine the finite-temperature behavior of 
the model. We can therefore speculate that the energy 
scales T+, T_ and A will reflect in analogous tempera- 
ture scales, leading to the scenario drawn in Fig. [1] and 
discussed in the Introduction. In this perspective, even 
though the critical region around T_ = is not stable 
at T = 0, and in fact in our calculation is replaced by 
superconductivity, it should become accessible by raising 
temperature T when T+ > T > Max{T_,Tc), denoted as 
NFL (non-Fermi-liquid) in the figure. If we follow a path 
in the T —U space that avoids the bell-shaped supercon- 
ducting region, we must find a crossover from a correlated 
Fermi-liquid phase to a pseudogap state passing through 
the NFL region. 

Finally, the region above T+ has a simple interpreta- 
tion in the impurity model, where it corresponds to the 
local moment regime above the Kondo temperature. In 
the lattice model, it presumably translates into a phase 
with very poor lattice-coherence, although its precise 
properties are difficult to foresee in the absence of ac- 
tual finite-T calculations. 



VI. MODELING THE DYNAMICS IN THE 
SUPERCONDUCTING PHASE 



An important result for the impurity model is that, 
as soon as one of the relevant symmetry breaking per- 
turbation is introduced, like, e.g., the hybridization 
suddenly the non- Fermi liquid behavior of the self-energy 
is replaced by a standard Fermi liquid one. In this case 
the low-frequency self-energy at the critical point changes 
from an imaginary constant to a conventional linearly 
vanishing function.—. We can therefore expect a similar 
regularization to occur as soon as superconductivity is 
allowed in the lattice model. This expectation is closely 
reminiscent of the mechanism taking place when s-wave 
superconductivity establishes in a disordered metal. Here 
scattering off impurities makes the quasiparticle lifetime 
1/r finite, which means that the normal state self-energy 
is finite and imaginary at zero frequency, just like in 
the above discussed non-Fermi liquid phase. Following 
Abrikosov, Gor'kov and Dzyaloshinskii,'*!., we write the 



self-energy as a 2x2 matrix whose diagonal entry En is 
the normal component, and the off-diagonal S12 is the 
anomalous (superconducting) contribution. We define 
the function ri{iujn) in the normal phase according to 

1 

iujn ~ Eii(icj„) = iujn + TT- sign(w„) = iojn ri{iujn)- (27) 

The onset of superconductivity regularizes this normal- 
state anomaly giving rise to a normal self-energy Sn lin- 
ear below a low-energy scale A and, at the same time, the 
anomalous self-energy E12 gets strongly enhanced with 
respect to a clean superconductor, namely^ 

zw„ - Ell {i^n) = itOn V + A2^ (28) 

El2(^c^„) = A77(V^^ + , (29) 

rj being the same function defined in Eq. (|27p . The 
two above equations express an important physical prop- 
erty, namely that non-magnetic disorder is a non pair- 
breaking perturbation in a BCS superconductor. Indeed 
Eqs. ([28|) and ([29|) imply a perfect cancellation in the s- 
wave Cooper channel between the wave-function renor- 
malization (the self-energy) and the vertex corrections 
brought by the impurities. This leads to the well-known 
"Anderson's theorem" stating that the value of Tc is 
independent of the concentration of non magnetic im- 
purities, provided the latter is low. This result follows 
immediately from the BCS gap-equation in the presence 
of an attractive coupling A^^ 



V 



A2 



(30) 



Tc is determined by solving ([30]) with A = 0, namely 
r]{iu!n) = 1-1- 1/(2t|w„|). One readily realizes that, sum- 
ming over momentum first, rj{iujn) disappears from the 
equation in the infinite bandwidth limit, leading to the 
same logarithmically singular sum over Matsubara fre- 
quencies as in the absence of disorder. 

We discussed previously that the emergence of SCS 
can be explained assuming that vertex corrections com- 
pensate exactly for the strong wave- function renormaliza- 
tion. It is then tempting to further pursue the analogy 
with dirty s-wave superconductors. Namely, we define 



lUJn 



Eii(icj„) 



(31) 



where Eii(ia;„) is the self-energy of the metastable nor- 
mal solution, using instead of •q{iujn), the more con- 
ventional notation Ziiujn) for the frequency dependent 
wave-function renormalization. We then assume that 
allowing for superconductivity leads to diagonal, En, 
and off-diagonal, E12, self-energy matrix elements, in the 
Nambu-spinor notation, given as in P5|) and (P^)) by 



iw„ - Ell («w„) 



A2 



(32) 
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£l2 {iuJn) = 



(33) 



13)1 . In Refs. I35ll36l the impurity model 



n I J 
2 



where A is the cut-ofF scale introduced by supercon- 
ductivity. It follows immediately that the introduc- 
tion of A restores a conventional Fermi-liquid behavior 
Sii —iuJn, not only when T_ = and Z{iujn) ^ \lO' 
but also in the pseudogap regime where Z{iujn) ~ \lU'_^^ 
Within this picture superconductivity regularizes the 
low-frequency behavior of the normal state self-energy 
both at the critical point, which is actually avoided by 
the onset of symmetry breaking, and in pseudogapped 
phase whose singularity is cut-off by A. As a consequence 
of such regularization, the low-energy pseudogap in the 
spectral function should change within the superconduct- 
ing phase into rather sharp quasi-particle peaks at the 
edge of the gap which finally disappears approaching the 
Mott insulator. 
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FIG. 8: (Color online) Fit of the impurity diagonal and ofF- 
diagonal self-energy, in the Nambu-spinor space, with the ex- 
pressions Eqs. p2|) and (|33p . 



Unfortunately a thorough comparison of numerical 
DMFT data with ^ and ^ is a very hard task, be- 
cause the phenomena we want to observe involve the ex- 
tremely small-frequency range (where the Fermi-liquid 
behavior would be recovered), which is really hard to 
study with present impurity solvers. For this reason, 
since the idea we propose is quite general, we have de- 
cided to verify the validity of our ansatz for the impurity 
model without self-consistency, postponing some selected 
DMFT comparison to the end of the section. 

An important difference is that no spontaneous sym- 
metry breaking is possible for the impurity model with- 
out self-consistency. Therefore, in order to check the 
validity of ([5^ and one has to add explicitly 

a symmetry-breaking perturbation to the Hamiltonian 



k.acr k.acr 



k a— 1 (J 



(34) 



was considered with U = 8, t± — 0.05 and different F's, 
in units of half the conduction bandwidth. In this model 
the role of the hybridization t±_ is twofold. On one hand 
it generates an exchange J± = At\/U able to drive the 
model through the critical point. At the same time, ij_ 
breaks the 0{2) orbital symmetry, turning the quantum 
phase transition into a crossover which has been shown to 
be quite sharp i^i^^ We note that, at particle-hole sym- 
metry, the two channels (jl]) and (fT8|) are perfectly equiv- 
alent, being related by the particle- hole transformation 



^2] 



'^2i ~^ "^21 (35) 
'^k2T ~' ~''k.2i '^k2i ~^ '^k,2T (3^) 

(37) 

where k and are particle-hole partners, ek = — ek, 
and V\^^ — V^. Therefore we can simply borrow the 
NRG data of Refs. I35ll36l and adapt them to our case of a 
superconducting symmetry-breaking term. In Fig. [5] we 
show the NRG data against our fit using Eqs. and 
with the normal self-energy that follows from the 



model DOS ([TBI). The fitting parameters T+, T_ and A 
are shown as function of F in Fig. [9l^^ We take the va- 




FIG. 9: (Color online) Fitting parameters T+, T_ and A of 
the impurity-model self-energy with the ansatz Eqs. (|32p and 
(|33p as a function of the hybridization F. 



lidity of Eqs. and ([55)1 for the impurity model as a 
strong support of their validity also in the lattice model 
even deep inside the pseudogap phase where, as previ- 
ously mentioned, our numerical data have not enough 
precision to reveal the very low-frequency structure. An 
interesting feature of the off-diagonal self-energy in the 
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impurity model is the strong frequency dependence, es- 
pecially in the pseudogap region where S12 is extremely 
peaked at very low frequencies. This is true also in the 
actual DMFT calculation. In Fig.[Tn]we plot, for different 
U's around the top of the SCS region, U ~ 0.82 W, the 
frequency dependent superconducting gap, defined by 



quasiparticle anomalous Green's function 



(38) 



whose zero-frequency extrapolation is the gap shown 
in Fig. [3] We note the rapid rise of A{iujn) below a 
frequency of order T^, confirming our prediction that, 
around its maximum, the gap is controlled by a single 
energy scale, T+. Therefore, even though the attraction 
J is instantaneous, strong retardation effects develop to 
avoid the large repulsion U and stabilize superconduc- 
tivity. It is just this strong frequency dependence of the 
superconducting gap ([55)1 that marks the difference be- 
tween a conventional Bose-Einstein condensation of pre- 
formed local pairs and our SCS phase, where pairs are 
highly non-local in time, see Fig. [TOl We also note that 
a characteristic energy scale shows up in the dynamical 
gap function even if our model does not involve any ex- 
change of bosons with a typical energy scale (at least not 
in an obvious way), as it has been instead proposed for 
the two-dimensional Hubbard model^. 
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FIG. 10: (Color online) Frequency-dependent superconduct- 
ing gap for several values of U. 

Finally, we conclude this section by discussing the pos- 
sibility to write down a gap-equation like Eq. ([30]) for 
our model in the SCS phase, which remains an open and 
suggestive issue. Starting from the Fermi-liquid side, our 
DMFT results seem to suggest that pairing essentially in- 
volves only the strongly renormalized quasiparticles. In 
this case a natural candidate for a gap equation would 
be 



(39) 



where is the renormalized scattering amplitude in the 
Cooper channel defined through Eq. (|12p while Fqp is the 



Fqp = ^F{iujn,ei^) 



(40) 



A first analysis shows that Eq. ([55]) correctly reproduces 
the order of magnitude of A also in the strongly corre- 
lated regime. This encourages us to use it to estimate 
the effective attraction A* such that Eq. ll39|) gives the 
actual DMFT value of A. We use ^ and (gD]) also in 
the pseudogap state, where the Fermi-liquid ansatz for 
A^,, Eq. (fT2|) . is not valid. Here Z is taken to be the 
low-energy spectral weight within a window of the order 
4J around the Fermi level. A^ that we obtain this way, 
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FIG. 11: Attractive coupling constant A, extracted through 
Eq. pO|l as the value that gives the correct DMFT value of 
A. 

shown in Fig. [11] is quite smooth as a function of U and 
it stays relatively small, in the range between 0.01 and 
0.02, up to the Mott transition. This shows that, even in 
the pseudogap regime, Eq. ([M[) can be satisfied without 
requiring a big coupling constant. 



VII. DISCUSSION AND CONCLUSIONS 

It is common wisdom that the emergence of local mo- 
ments out of the incipient Mott localization must nec- 
essarily bring to an enhanced magnetic response and 
eventually to a magnetic order that may appear already 
in the metallic phase adjacent the Mott insulator. It 
is equally conceivable that, under these circumstances, 
a system may become easily unstable to lattice distor- 
tions. Moreover, on the brink of charge localization, also 
the orbital momentum, quenched by the hopping deep 
in the metallic phase, may re-emerge making spin-orbit 
coupling effective as if atoms were isolated. This is the 
typical phenomenology of magnetic Mott insulators ex- 
emplified by the prototype Mott-Hubbard system, the 
Y^Q^ iSAe These properties can be interpreted within 
a Landau Fermi-Liquid framework only by invoking an 
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almost perfect cancellation between the large wave func- 
tion renormalization associated to the Mott transition 
and the vertex corrections in peculiar channels. Physical 
intuition suggests that such channels should be primarily 
those acting on degrees of freedom orthogonal to charge, 
like spin and orbital momentum. 

However, while the above arguments are quite nat- 
ural, if not obvious, for particle-hole instabilities like 
magnetism, they sound much less trivial in connection 
with superconducting particle-particle instabilities. The 
main difference is that superconductivity implies phase- 
coherence, which can not survive charge localization. 
Nevertheless, while phase-coherence requires pairing, the 
opposite is not true, since a Mott insulator can be formed 
by incoherent pairs. An example is a valence-bond crys- 
tal formed by an ordered array of tightly bound singlets 
sitting on nearest-neighbor bonds, and, seemingly, also 
a Mott insulator made of resonating valence bonds - 
the RVB scenario originally proposed by Anderson^^ for 
high-Tc superconductors. 

Another possibility is to form local singlets exploiting 
the orbital degrees of freedom, that play a role similar 
to the bonds of the previous examples. This is realized 
in our model ((T|), and in model for fuUerenes. In the 
fuUerenc family, tetravalent alkali-doped Cgo may be re- 
garded as the parent Mott insulating compound^^'^^ that 
turns superconducting upon doping, the trivalent mate- 
rials. The former compounds are indeed non-magnetic 
Mott insulators where the four electrons occupying the 
LUMO of each molecule bind into a non degenerate spin- 
singlet configuration because of Jahn- Teller effect.— 

In all these examples, the pairing implicit in the Mott 
state is perfectly compatible with strong repulsion, hence 
it is not surprising that, moving away from the Mott insu- 
lator, superconductivity appears. The general conditions 
are elucidated by our analysis: pairing must correlate de- 
grees of freedom orthogonal to charge. For instance, in 
the t-J model for cuprates pairing is provided by the 
spin superexchange, that is unaffected by the strong re- 
pulsion constraint of no-double occupancy, while it com- 
petes with the hopping. The former favors configurations 
in which two singly-occupied nearest neighbor sites are 
bound into a singlet state, while the latter prefers a demo- 
cratic occupancy, in which the singlet is equally probable 
as any of the triplet states. Also in our two-orbital model 
([T|) close to the Mott transition, J competes with the hop- 
ping rather than U . At half-filling the hopping favors an 
equal occupation of all the states with two electrons per 
site, while J breaks this degeneracy in favor of the singlet 
configuration ([5]). Right the same competition emerges 
in the impurity model (fO)) between the Kondo tempera- 
ture and J. The reason why J <C i eventually prevails in 
all the above examples is that the hopping suffers from a 
very severe wavefunction renormalization close to a Mott 
transition, while J does not. 

From these observations it appears evident that the 
physics of our model ((T]) and that of the t-J model for 
cuprates within the RVB scenarici^ii^'iSi'^ share com- 



mon features; our inverted exchange playing on-site a 
similar role as the nearest-neighbor exchange in the t-J 
model. Quite obviously J being on-site or on a link is ex- 
pected to introduce relevant differences as far as the mo- 
mentum structure is concerned (the most evident being 
s-wave pairing versus d-wave pairing) . There are however 
other relevant differences, at least regarding the interpre- 
tation of the various phases. In the RVB scenario uncov- 
ered by slave-bosoni^ii^iSSiM and variational^^i^i^i^iS 
approaches, a lot of emphasis is placed on the pseudo- 
gap regime close to the Mott transition to explain super- 
conductivity.^^ In our model superconductivity is instead 
the low-temperature response to the instability in the 
crossover region between the pseudogap and the Fermi- 
liquid. Another important difference is that our pseu- 
dogap state is a stable thermodynamic phase with zero 
entropy at zero temperature, while in the RVB scenario 
it appears as a fluctuating mixture of competing orders, 
unable to survive down to zero temperature. 

Our results look also different from QCP theories for 
the cuprates^iSiiiS. Within the QCP approach super- 
conductivity arises because of the critical fluctuations 
around a true quantum critical point that separates sta- 
ble zero-temperature phases. This critical point may 
even be inaccessible in the physical space of parame- 
ters, and exist only in an hypothetical enlarged space, 
yet it is believed to influence the physical system not 
only at finite but also down to zero temperature. In con- 
trast, in model ([T]) the normal Fermi-liquid metal, the 
pseudogap phase, the superconductor and other possible 
symmetry-broken phases are all equally legitimate out- 
comes of an underlying competition that reveals itself 
only in the high-temperature non-Fermi-liquid crossover 
region, NFL in Fig. [T] From this point of view, the dif- 
ference between the phase diagram Fig. [7] and the QCP 
scenario is the same emphasized recently by Andersoni^ 
An obvious and pertinent criticism to our claim is that, 
upon applying a sufficiently high magnetic field able to 
completely wash out superconductivity, the unstable crit- 
ical region should move down to zero temperature and 
transform into a true quantum critical point thus recov- 
ering the QCP scenario. We note however that, at its 
maximum, the superconducting gap is only controlled by 
the energy scale T+, which is also the bandwidth of the 
low-energy incoherent excitations. For this reason we 
believe, although we can not prove it, that a magnetic 
field so strong to completely suppress superconductivity 
would drive the model into a phase with very poor lattice 
coherence, rather then revealing the critical point. 

We believe that our results provide a new perspec- 
tive about the relation between a pseudogap and su- 
perconductivity, which crucially differs from previous ap- 
proaches. Rather than being in competition, the pseu- 
dogap and the superconducting gap turn out to be much 
more "compatible" one another than in most theoreti- 
cal approaches that we are aware of. A central aspect, 
discussed in Sec. I VII is that the large normal self-energy 
responsible of the pseudogap state is regularized by the 
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onset of superconductivity, in contrast with approaches 
in which the same (highly anomalous) normal self-energy 
is assumed in the normal and superconducting states^. 
Moreover, we have shown that in model ^ there is no 
need of a large coupling constant to overcome the lack 
of low-energy density of states in the normal phase, as 
it happens instead in approaches where pairing takes 
place starting from a normal state with mean-field like 
pseudogapsi^iS. 

Another important outcome of our calculation is the 
natural appearance of two energy gaps with different be- 
havior as a function of the distance from the Mott in- 
sulator (see Fig. [1]): the pseudogap scale T_ increases 
approaching the Mott state, where it is largest, while the 
superconducting coherence scale vanishes as the Mott in- 
sulator is approached. This is clearly reminiscent of the 
two energy scales (gaps) observed in the cuprate a- -' ' - "" . 

Finally, we want to emphasize that the physics we have 
unraveled seems to be more general than the specific toy- 
model ([H). Indeed, we have shown that the emergence 
of a reinforced superconductivity upon approaching the 
Mott transition can be explained in the model ([T]) by 
simple Fermi-liquid arguments that do not depend on 
the specific model. They simply state that any scatter- 
ing channel involving degrees of freedom orthogonal to 
charge should strengthen near a Mott transition, irre- 
spectively whether being particle-particle or a particle- 
hole channels. This is for instance obvious in one dimen- 
sion because of the dynamical separation of the charge 
from the other degrees of freedom, and explains why the 



phase diagram of U) in one dimension does not differ 
qualitatively from Fig. [7p^ In addition, the impurity- 
model analogy, that we have exploited in this work, indi- 
cates that some kind of pseudogap phase - bulk counter- 
part of an impurity unscreened regime - might be ubiq- 
uitous near a Mott transition, especially if this transition 
is continuous. Indeed it is obvious that a metal should 
suppress spectral weight around the chemical potential 
to smoothly connect to a realistic zero-entropy insulator, 
and, in the absence of symmetry breaking, we do not see 
any other way but opening a pseudogap. This would in 
turn imply also the existence of a high-temperature cross- 
over regime that separates the Fermi-liquid phase from 
the pseudogap one and actually contains the seeds of the 
symmetry broken phase that eventually emerge at low 
temperature.'^^ However, it often happens that the Mott 
transition is first order, in particular when accompanied 
by a lattice distortion, like in ¥203.^^ In this case both 
the pseudogap and the critical cross-over regions might 
become unaccessible and a discontinuous transition occur 
directly from the Fermi-liquid metal to the Mott insula- 
tor. 



Acknowledgments 

We are grateful to Lorenzo Dc Leo for providing us 
with his numerical renormalization group data on the 
impurity model. 



^ C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, 

Phys. Rev. Lett. 71, 3202 (1993). 
2 R. Haussmann, Z. Phys. B 91, 291 (1999). 
^ N. Andrenacci, A. Perali, P. Pieri, and G. C. Strinati, Phys. 

Rev. B 60, 12410 (1999). 
* B. Kyung, S. Allen, and A.-M. S. Tremblay, Phys. Rev. B 

64, 075116 (2001). 
^ A. M. Tsvelik and A. V. Chubukov, Phys. Rev. Lett.) 98, 

237001 (2007), and references therein. 
^ Y. Wang, L. Li, and N. P. Ong, Phys. Rev. B 73, 024510 

(2006). 

A. Kanigel et al., Nat. Phys. 2, 447 (2006). 
® C. M. Varma, Phys. Rev. B 55, 14554 (1997). 
^ C. Castellani, C. Di Castro, and M. Grilli, Z. Phys. B 103, 

137 (1997). 

1° J. L. Tallon and J. W. Loram, Physica C 349, 53 (2001). 
" A. J. Millis, Phys. Rev. B 48, 7183 (1993). 

P. Coleman, C. Pepin, Q. Si, and R. Ramazashvili, J. 

Phys.: Condens. Matter 13, 723 (2001). 

Q. Si, Adv. in Solid State Phys. 44, 253 (2004). 

P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 

78, 17 (2006), and references therein. 
^5 S.-C. Zhang, Science 275, 1089 (1997). 

P. W. Anderson, Science 235, 1196 (1987). 

A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 

Rev. Mod. Phys. 68, 13 (1996). 

M. Capone, M. Fabrizio, C. Castellani, and E. Tosatti, 



Phys. Rev. Lett. 93, 047001 (2004). 
1^ G. Kothar and J. Liu, Phys. Rev. B 38, 5142 (1988). 

Y. Suzumura, Y. Hasegawa, and H. Fukuyama, J. Phys. 

Soc. Jpn. 57, 2768 (1988). 
2^ O. Gunnarsson, Rev. Mod. Phys. 69, 575 (1997). 

M. Capone, M. Fabrizio, and E. Tosatti, Phys. Rev. Lett. 

86, 5361 (2001). 

M. Capone, M. Fabrizio, C. Castellani, and E. Tosatti, 
Science 296, 2364 (2002). 

E. Plekhanov, S. Sorella, and M. Fabrizio, Phys. Rev. Lett. 
90, 187004 (2003). 

R. Assaraf, P. Azaria, E. Boulat, M. Caffarel, and P. Le- 
cheminant, Phys. Rev. Lett. 93, 016407 (2004). 
'^^ F. F. Assaad, Phys. Rev. B 71, 075103 (2005). 

A. Paramekanti and J. B. Marston, J. Phys. Cond. Matt. 
19, 125215 (2007). 

M. Fabrizio and E. Tosatti, Phys. Rev. Lett. 94, 106403 
(2005). 

^8 M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 
(1994). 

^° M. Capone, M. CivelU, S. S. Kancharla, C. Castellani, and 
G. Kothar, Phys. Rev. B 69, 195105 (2004). 
W. Metzner and D. VoUhardt, Phys. Rev. Lett. 59, 121 
(1987). 

'^^ R. Micnas, J. Ranninger, and S. Robaszkiewicz, Rev. Mod. 
Phys. 62, 113 (1990). 

A. Toschi, M. Capone, and C. Castellani, Phys. Rev. B 72, 



14 



235118 (2005). 

M. Fabrizio, A. F. Ho, L. Do Loo, and G. E. Santoro, Phys. 
Rev. Lott. 91, 246402 (2003). 

L. De Leo and M. Fabrizio, Phys. Rev. B 69, 245114 
(2004). 

M. Ferrero, L. De Leo, P. Lecheminant, and M. Fabrizio, 
J. Phys. Cond. Matt. 19, 433201 (2007). 
H. R. Krishna-rnurthy, J. W. Wilkins, and K. G. Wilson, 
Phys. Rov. B 21, 1003 (1980). 

H. R. Krishna-rnurthy, J. W. Wilkins, and K. G. Wilson, 
Phys. Rev. B 21, 1044 (1980). 

B. A. Jones and C. M. Varma, Phys. Rev. B 40, 324 (1989). 
*° 1. Affleck, A. W. W. Ludwig, and B. A. Jones, Phys. Rev. 
B 52, 9528 (1995). 

A. A. Abrikosov, L. P. Gor'kov, and I. E. Dzyaloshinski, 

Methods of Quantum Field Theory in Statistical Physics 

(Dover, New York, 1963). 
*2 P. W. Anderson, J. Phys. and Chem. Sohds 11, 26 (1959). 
''^ Since the impurity model contains explicitly the symmetry 

breaking term, the solf-onergy that wc fit is actually the 

off-diagonal component of the inverse Green's function. 

A. Macridiu, M. .Jarroll, T. Maior, and D. J. Scalapino, 

Phys. Rev. Lett. 99, 237001 (2007). 

W. Bao, C. Broholm, G. Aeppli, S. A. Carter, P. Dai, T. F. 
Rosenbaum, J. M. Honig, P. Metcalf, and S. F. Trevino, 
Phys. Rev. B 58, 12727 (1998), and references therein. 
L. Paolasini, C. Vcttior, F. do Bergevin, F. Yakhou, 
D. Mannix, A. Stunault, W. Neubeck, M. Altarelh, 



M. Fabrizio, P. A. Metcalf, et al., Phys. Rev. Lett. 82, 

4719 (1999). 

M. Fabrizio and E. Tosatti, Phys. Rev. B 55, 13465 (1997). 
G. Baskaran, Z. Zou, and P. W. Anderson, Sol.State 
Comm. 63 (1987). 

C. Gros, Phys. Rev. B 38, 931 (1988). 

^° S. Sorella, G. B. Martins, F. Becca, C. Gazza, L. Capriotti, 
A. Parola, and E. Dagotto, Phys. Rev. Lett. 88, 117002 

(2002) . 

®^ A. Paramckanti, M. Randeria, and N. Trivedi, Phys. Rev. 
Lett. 87, 217002 (2001). 

P. W. Anderson, P. A. Lee, M. Randeria, T. M. Rice, 
N. Trivedi, and F. C. Zhang, J. Phys. Condens. Matter 
16, R755 (2004). 

D. A. Ivanov and P. A. Lee, Phys. Rev. B 68, 132501 

(2003) . 

P. W. Anderson, Physica B 318, 28 (2002). 
55 K.-Y. Yang, T. M. Rice, and F.-C. Zhang, Phys. Rev. B 

73, 174501 (2006). 
5« P. Nozicrcs and F. Pistolosi, Eur. Phys. J. B 10, 649 (1999). 
5'^ L. Bonfatto, S. Caprara, and C. D. Castro, Eur. Phys. J. 

B 17, 95 (2000). 
5® M. Le Tacon, A. Sacuto, A. Georges, G. Kothar, Y. Gallais, 

D. Colson, and A. Forget, Nat. Phys. 2, 537 (2006). 
5® K. Tanalca, W. Lee, D. Lu, A. Fujimori, T. Fujii, Risdiana, 

I. Terasaki, D. Scalapino, T. Devereaux, Z. Hussain, et al., 

Science 314, 1910 (2006). 



